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Abstract 

' With the aid of concrete examples, we consider the question of whether, in the presence 

of conformal curvature, a conformal geodesic can become trapped in smaller and smaller 
sets, or phrased informally: are spirals possible? We do not arrive at a definitive answer, 
but we are able to find situations where this behaviour is ruled out, including a reduction of 
the conformal-geodesic equations to quadratures in a specific lum-confonnallv Hat iiieine. 
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1 Introduction 



Conformal geodesies are a family of curves naturally defined in a conformal manifold of any 
signature. They play a role with respect to the conformal structure similar to that played by 
the geodesies of a metric, and the definition has been around at least since the 1950s, [T], but 
much less is known about them than about metric geodesies. For example there are unanswered 
| questions related to completeness of the form: can a conformal geodesic enter and not leave 

smaller and smaller sets? A specific version of this question due to Helmut Friedrich [2] is: if a 
conformal geodesic T enters every neighbourhood of a point p, does T necessarily pass through p 
(with a finite limiting velocity and acceleration)? There are similar and related open questions in 
Schmidt's construction of the c-boundary [3]. A metric geodesic cannot be trapped in smaller and 
smaller sets and a proof can be based on the existence of geodesically-convex neighbourhoods: 
if a geodesic enters a geodesically-convex ball then it must leave it too. Without a similar 
notion of convexity for conformal geodesies, it is not clear what the answer to the corresponding 
question should be. In conformally-flat metrics, such as Euclidean or Minkowski spaces, there 
are conformal geodesies which are closed plane circles of arbitrarily small radius. It is usually 
thought that these will fail to close up if some conformal curvature is added (and we shall see 
examples confirming this below), but whether the perturbed conformal geodesies can then spiral 
in in some fashion, either to an arbitrarily small circle or to a fixed point, is not known. 

In this article we approach the question of spirals with the aid of some examples, including 
a reduction to quadratures of the conformal geodesic equation in the Riemannian or Lorentzian 
metrics on Nil, one of Thurston's geometries (see e.g. [4]; I believe this is the first published 
instance of a reduction to quadratures of the conformal geodesic equation in a non-conformally 
flat metric but, possibly just because this case is integrable, there are no spirals). We do not 
solve the problem of the existence of spirals in general but we are able to establish nonexistence 
in certain cases, and we find various other properties of conformal geodesies. 

The object of study then is a conformal manifold (M, [g]), where [g] is a conformal equivalence 
class of nondegenerate, say smooth, metrics which can have any signature, though we shall 
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usually restrict to Riemannian or Lorentzian (where this is taken to mean one negative eigenvalue 
for the metric). 

Make a choice of metric g ao in the conformal class, with corresponding metric covariant 
derivative V Q . Then a conformal geodesic is a curve T with tangent vector v a and a one-form 
b a given along T and satisfying the system: 

v c V c v a = -2(v c b c )v a + g ac b c (g ef v e vf) (1) 
v c V c b a = (v c b c )b a -^g ac v c (g^b e b f )+L ac v c (2) 

where, in dimension n, L a b = ( n l 2 ) (R<tb ~ 2(n-i) ^ ab )' * s ^ ne ^ cc i tensor and i? is the 
scalar curvature. Other equivalent definitions of conformal geodesic are possible (see [5]). In 
indefinite signatures, the definition still makes sense for null v a when ([1]) reduces to the null 
geodesic equation: null geodesies are conformal geodesies, and the only null conformal geodesies 
are null geodesies. Henceforth we shall assume that the conformal geodesies considered are not 
null. 

The conformal geodesic equations are conformally-invariant in that, under conformal rescal- 

ing 

9ab -> 9ab = ^9abi (3) 

with the usual corresponding change to V a , they take the same form with 

v a = v a (4) 

b a = b a ~ Y a (5) 

where, as usual, T a = f2 _1 V a fi. 

A conformal geodesic admits a preferred parameter r defined up to choice of origin by 

v c V c t = 1. (6) 

There is a reparametrisation freedom of Mobius transformations in r: 

ar + b 



ct + d' 



(7) 



so that r may be called a projective parameter, which must be accompanied by a transformation 
of v a and b a : 

v a = {cT + dfv a (8) 
k = b a + fg ac v c (9) 

where 

f = -2(g ef v e vfy 1 c(cT + d)- 1 . 

This transformation draws attention to a familiar problem in the study of conformal geodesies: 
it is possible for v a to vanish and b a to be singular at a regular point both of the manifold and 
of the curve, by choice of projective parameter. One way to avoid this problem is to introduce 
a third-order form of the equations. How this is done depends slightly on the signature. In the 
Riemannian case or the Lorentzian case with a space-like velocity, introduce the unit tangent 
u a = x _2 ^ a , where x 4 = 9efV e v^ > 0, then ([T]) can be solved for b c in terms of the acceleration 
a c = u a V a u c as 

b c = a c + (u a b a )u c , (10) 
where u a b a = —2x~ 1 u a V a X- Use this in ([2]) to find 

u c V c a b = u b (-g ef a e a f - L ef u e u f ) + L b c u c , (11) 
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which is the third-order form of the conformal geodesic equation, and 

u c V c (u a b a ) = l -{g eJ a e af + (u% a ) 2 ) + L ef u e u^ 

which can be written as 

X = -\{9e f a e af + 2L ef u e uf) X , (12) 

where the overdot denotes differentiation with respect to proper length. We shall always call 
the parameter t, even when it is proper length. From ([6]) and the definition of u a we have 

d T = v a d a = X 2 u a d a = x 2 d u 

so that dr = dt/x 2 ■ 

Having solved (|lip . to recover a representative (v a , b a ) we solve (|12p for a choice of x when 
v a = x 2u ° \ an d b a is obtained from fjlOf) . The different choices of x correspond to the different 
choices of the projective parameter r. For example, if we choose two solutions XI1X2 of (|12p 
with unit Wronskian: 

X1X2 ~ X2X1 = 1, so that ^- (—) = \ and X\ 2 dt = d ( — J , 

dt \XiJ Xi \XiJ 

so that an allowed choice for r is Xilxx- Thus zeroes and poles of the projective parameter are 
tied to zeroes and poles of the solutions of (fT2"j) , 

For the time-like Lorentzian case, introduce unit tangent u a = x -2 ^ where now x 4 = 
—g e fV e vf, then (pQ) can be solved for b c in terms of the acceleration a c = u a V a u c as 

b c = -a c - (u a b a )u c , (13) 

where u a b a = —2x~ 1 u a V a x- Use this in ([2]) to find 

u c V c a b = u\g ef a e a f - L ef u e u f ) - L b c u c , (14) 

and 



X = \{9efa e a f -2L ef u e u f )x- (15) 



Now for existence we have: 



Theorem 1.1 A conformal geodesic, that is a solution of or J 1J$ , will always exist for 
some interval (— ti,ii) in t, given data at time t = consisting of an initial position, an initial 
unit velocity u a and an initial acceleration a a orthogonal to u a . Thus a conformal geodesic can 
always be continued in the Riemannian case as long as a a remains finite, and in the Lorentzian 
case as long as both u a and a a remain finite. 

Proof 

This follows at once from Picard's Theorem. 

□ 

It follows as a corollary that spiralling can only occur in the Riemannian case when the 
acceleration diverges or in the Lorentzian case when either the velocity or the acceleration 
diverges (though presumably one cannot happen without the other). For conformal geodesies 
in flat space or in an Einstein manifold, the magnitude of the acceleration is constant. Thus 
spiralling cannot happen in a Riemannian Einstein space, and it seems likely that this is also 



3 



true in Lorentzian Einstein spaces. In the presence of more general curvature, the acceleration 
need not be constant since we have, from (fl~Tj) or (fTl"|) . that 

^(g ef a e a f ) = ±2a e u f L ef 

with the sign depending on the signature. If we set a e = qw e where w e is a unit vector and 
q > 0, then 

q = ±w e n,f L e f, 

where the plus sign is for Riemannian or Lorentzian space-like and the minus sign is for 
Lorentzian time-like. Now for Riemannian, in a region in which the curvature is bounded, the 
magnitude q of the acceleration has bounded time-derivative and so the acceleration can only 
diverge after infinite t, that is after infinite proper distance, and so spiralling could only occur 
after infinite proper distance. Furthermore, if the magnitude of the acceleration does diverge 
then by a simple comparison argument applied to ()12p . any admissible x will have infinitely 
many zeroes, and so any choice of projective parameter r will pass infinitely often through zero 
and infinity. 

The Lorentzian cases are different, for example both the velocity and the acceleration can 
become unbounded in finite time, as we shall see from examples below. 

In the rest of this article, we briefly review in Section 2 the behaviour of conformal geodesies 
in three-dimensional Einstein spaces. Much of this carries over to conformally-flat spaces of any 
dimension and is quite untypical of the general case. 

In Section 3, we consider some very simple axisymmetric metrics in which it is possible to 
solve enough of the conformal geodesic equations to show that there are some which close up 
into circles but plenty more which do not: the widespread intuition that conformal curvature 
typically prevents circles from closing is vindicated. What makes the calculation tractable is 
that planes of constant z in these metrics are totally-geodesic for conformal geodesies (call this 
property totally- conformally geodesic). It is also possible to show in these examples that the 
conformal geodesies in these 2-dimensional totally-conformally geodesic surfaces cannot spiral 
into the origin or into arbitrarily small limit cycles round the origin - certain kinds of spiralling 
are ruled out. 

In Section 4, we reduce the conformal geodesic equations in a particular metric to quadra- 
tures. The metric is a Riemannian or Lorentzian version of Nil (which arises in Thurston's 
classification of geometries, see e.g. [4]). The reduction is possible because of the existence of 
two constants of the motion. One comes from a Killing vector with a particular algebraic rela- 
tion to the Ricci tensor; the other arises because the Ricci tensor is a Killing tensor. One of the 
constants serves to bound the norm-squared of the acceleration so that, in the Riemannian case, 
the acceleration cannot blow up and any conformal geodesic can be continued forever - there are 
no spirals in this example. In the Lorentzian case the indefiniteness of the signature means that 
a bound on the norm-squared of the acceleration does not give a bound on the components of 
acceleration and one can find explicit examples with the acceleration blowing up exponentially 
in time (which can already happen in Minkowski space) or in finite time (which cannot). Again, 
there are no spirals. 

In Appendix A we show that the conformal geodesic equations in the Berger sphere can also 
be reduced to quadratures. This is not too surprising given the result for Nil because the metrics 
have a number of similarities, including the Killing tensor property of the Ricci tensor. 

In Appendix B we write out the conformal geodesic equation in the Schwarzschild metric, 
a Lorentzian Einstein metric. The novelty here is that there are two constants of the motion 
arising from the Killing spinor, though none from the Killing vectors so that we do not have 
enough constants for integrability. 
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2 The cases of S 3 ,E 3 ,H*, M 3 



By this we mean the three, three-dimensional Riemannian Einstein metrics of positive, zero and 
negative curvature respectively, and three-dimensional Minkowski space. These cases are well 
understood but it is useful to review them for what follows. All are conformally-flat, so that 
the conformal geodesies in S 3 and H 3 can be obtained from those in E 3 , and all are Einstein, 
so that the term in L ab drops out of (fTTj) or (fT3|) . We're left with 

u c V c a b = -eu\g ef a e af), (16) 

where e = — 1 for Lorentzian time- like and e = 1 for Riemannian or Lorentzian space-like (so 
that e = g ab u a u b ). We see at once that g e fa e af is constant. Since this can be zero, we see that 
all metric geodesies are actually conformal geodesies. This will not be true with more general 
curvature: 

Remark: By inspection of (jlip or (I14p . we see that only a metric geodesic whose tangent vector 
is an eigenvector of the Ricci tensor will also be a conformal geodesic. 

To find more constants of the motion, suppose K a is a Killing vector, so that 
M ab := V a K b = V [a K b] , V a M bc = ^R{g ab K c - g ac K b ), 
where R is the Ricci scalar, and we've used the expression 

Rabcd = -j:R{9ac9bd ~ 9ad9bc), 
O 

for the Riemann tensor, which is valid in these metrics. Now consider the quantity 

Q := M ab u a a b - l -RK a u a . (17) 
b 

It is easy to see that this is constant, and this observation yields six constants of the motion, 
plus g e fa e a-f and g e fU e u^ for S 3 and H 3 . For E 3 and M 3 , where R = 0, only the 3 rotations give 
non-trivial constants like this (the translations give zero), and these correspond to the constant 
bivector U[ a a b ] which defines the plane of the orbit. 

To find the explicit solution for the conformal geodesies in E 3 , when e = 1, introduce 
A 2 = g e ja e af then (fT6j) becomes 

u a = -A 2 u a . 

If A = then the solution is a metric geodesic. If A ^ then this is solved by 

u a = UQCOs(At) +A~ 1 Qq sm(At), 

where = u a (0) and ag = a a (0); the position x a {t) follows by integration: 

x a = + A~ l ul sin(At) + A~ 2 al{l - cos(At)). (18) 

The conformal geodesies are planar circles and straight lines. 

To find the explicit solution in M 3 , in the time-like case (so e = — 1) a a is orthogonal to u a 
and so is space-like. Now g e ja e a^ = A 2 and (fTBT) is solved by 

u a = Uq cosh(At) + A~ l al sinh(At), 

whence again x a (t) follows by integration. This is the constant-acceleration world-line famil- 
iar from courses on special relativity: the scalars g ab u a u b and g ab a a a b are constant but the 
components of both velocity and acceleration are unbounded in a constant basis. 
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The Lorentzian space-like case can have a a time-like, space-like or null. The assumption of 
time-like leads to another hyperbola, and space-like gives another circle, so the only new case is 
null, when a a = a§ is a constant null vector and 

x a = xl + tug + -t 2 ag. 

This is a parabola lying in the null hyperplane g a bO-o( xb ~ x o) = 0- 

Solutions in 5 3 or H 3 can be obtained from those in E 3 by conformal rescaling. They turn 
out to be intersections of S s embedded in the standard way in E 4 or H s embedded in the 
standard way in M 4 by linear subspaces in the ambient space, so they are circles or hyperbolae. 

3 An axisymmetric example 

We can get enough information to throw some light on our questions from a simple axisymmetric 
example. First note: 

Theorem 3.1 Any trajectory of a (non-null) hypersurf ace- orthogonal Killing vector is a con- 
formal geodesic. 

Proof 

The condition of hypersurface-orthogonality means 

K [a V b K c] = 0, 

so that, since K a is a Killing vector, 

V a K b = V [a K b] =2V [a K b] , 

for some vector V a which can be assumed to be orthogonal to K a . 

Suppose 7 is a trajectory of the Killing vector and set K a K a = ell 2 where e = ±1, to allow 
for different signatures and U ^ by assumption. The unit tangent vector to 7 is u a = U~ 1 K a 
and the acceleration is 

a a = u h V b u a = -eV a , 
using what we know. Recall the Killing vector identity 

V a Vb.K" c = R bcac [K d , 

so that 

R cd K d = K c (V a V a + V a V a ) - K b V b V c . 
and it is easy to check that (jlip is satisfied. 

□ 

Remark: There's no reason to expect this result to hold for Killing vectors which are not 
hypersurface-orthogonal, and examples where it fails will be seen below, in Section 4.1. 

Now we'll do the example. Suppose the metric is 

ds 2 = dr 2 + F 2 (r)# 2 + dz 2 , (19) 
with a regular axis at r = 0, and consider conformal geodesies lying in z = 0. Introduce x by 

r = cosx, Fcp = smx 
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to take care of the normalisation of the velocity (since i = 0) , then the velocity is 



u = cos x e i + sin x&2 

in terms of orthonormal vectors 

ei = d r , e 2 = F~ 1 d (j} . 

We claim the acceleration is then 

a = q(— sin^ei + cosxe 2 ), 

where 

q = x + sinxF' /F, 

and the prime means d/dr. Because of the dimension of the surface z = 0, the Ricci term drops 
out of the conformal geodesic equation (|lip and we're back to (116p . The remaining equation is 
equivalent to the constancy of g a i,a a a b or equivalently the constancy of q, which may be written 

q = x + tan xF ' /F = constant, 

where we've used f = cosx- Introduce 

,. r 

G(r) = / F(s)ds, 
Jo 



then this integrates as 
from which we obtain 



Fsinx — qG = C = constant, 

. 2 (C + qG) 2 

r =1 ^ =-f( r )- 



For a regular axis at r = we need F/r — > 1 as r — > 0, so also 2G/r 2 — > 1. Thus if q and 
C are both nonzero then f(r) is negative for small r, and there is a lower limit for r on the 
conformal geodesic, say T. Behaviour at large r depends on what conditions we impose there. 
For asymptotic flatness we would want the same limits at large r as at small r and then, if q 
and C are both nonzero, f(r) is negative for large r, so that T is confined between finite and 
nonzero upper and lower limits in r. We also have the orbit equation as 

C + qG 

-dr, 



F(F 2 - (C + qG) 2 ) 1 / 2 

so that the conformal geodesies do not close in general, but rather fill out an annulus in the 
plane. This is a confirmation of the common belief noted in the Introduction. 

The special case of q = is actually a metric geodesic, which can escape to infinity in r. The 
special case of C = is a conformal geodesic passing through the origin. If q ^ 0, it will have 
an upper limit in radius r to which it will return after each pass through the origin - this is not 
a spiral but something more like a rosette. Indeed careful choices will give conformal geodesies 
which are closed rosettes, like the curve r = asin2(p. 

From Theorem 13. 1\ it must be possible to choose q and C so that any fixed value of r, say 
r = a, is a conformal geodesic. This requires f(a) = = f'(a) which in turn requires 

q = ±F'(a)/F(a), C = + {{F{a)f - F'(a)G(a)) /F(a). 

Now one might wonder whether there could be conformal geodesies which spiralled in asymp- 
totically in time to these circles. The existence of such conformal geodesies for arbitrarily small 
a would provide a negative answer to the question of Friedrich mentioned in the Introduction. 
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Remark: While this may be possible for some a, it cannot be arranged for arbitrarily small a. 
To see this calculate 

f"(a) = 2 (F(a)F"(a) - (F'(a)f) /{F{a))\ 

at a circular orbit. At small enough a this must be negative since F — > 0, F' — > 1 as a — > 0. 
Thus small circular conformal geodesies are stable, but to be the limit cycle of a spiral, a small 
circle would need to be unstable (or at worst marginally stable). 

Spirals can also be ruled out in a slightly more general case: replace (fTUj) by 

ds 2 = dr 2 + F 2 (r, z)d(j) 2 + H 2 (r, z)dz 2 , (20) 

but where we'll suppose F and H are even in z so that the plane z = is totally-geodesic for 
metric geodesies. We shall see below that conformal geodesies can also be confined to this plane. 
The velocity is 

a = f, f3 = F<fi, so that a 2 + /3 2 = 1. 
Then the acceleration is given by 

F' ■ F' 

ai = a-0 2 — , a 2 = $ + a(3—, 

retaining prime for d r , and the conformal geodesic equations turn out to be 

F 1 9 9 

di = Pci2— — \a\ a + a/3 Q, 
F 

pi 

«2 = —Pa\— - \a\ 2 (3 - a 2 f3Q, 

r 

where 

_FT_ F\^_ F'H' 

Q ~~1T + f1p + Th' 

which is a combination of curvature components, and in particular must be bounded as r — > 0. 
Since the velocity is a unit vector we must have 

a\a + (Z2/3 = so a\ = qf3, a 2 = —qa, 

for some q, when in fact q 2 = \a\ 2 . Now the conformal geodesic equations collapse to 

q = a/3Q = f(3Q, 

so dq/dr = (3Q which is bounded for small r. Therefore one cannot have a conformal geodesic 
spiralling in to the origin since along such a curve the magnitude of the acceleration must be 
unbounded, while this equation would force q to be bounded. 

It is also possible to check in this more complicated example that while the axial Killing 
vector trajectories provide small circular conformal geodesies, the smallest circles are again 
stable and cannot be the future limit cycles of spirals. 

We've just seen that the plane z = is totally-conformally-geodesic in that conformal 
geodesies could be confined to it. It is worth noting the following: 

Theorem 3.2 The necessary and sufficient condition for a hypersurface E to be totally-conformally- 
geodesic is the vanishing of the second- fundamental form of S . 
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Proof 

Let the unit normal to S be N a , then the covariant derivative of N a takes the form 

V a N b = N a A b + K ab , 

where 

N a A a = = K ab N b , K [ab] =0, 
so K ab is the second fundamental form of E. Start a conformal geodesic in S so that initially 

qi ■= N a u a = 0, q 2 := iV a a a = 0. 

Calculate 

qi=q 2 + K ab u a u b + gi(n a A a ), 

and 

42 = -qi(9a b a a a b ) + 7V a L afeU fe + u a a b (K ab + N a A b ). 

Now uniqueness of solution of this system will force qi = = q 2 for every initial choice provided 
K ab u a u b = = N a L ab u b for every allowed u a . This certainly requires K ab = 0. Suppose that 
holds and recall the Gauss-Codazzi equation: 

R bc h c a N b = D b K b - D a K, 

where h ab = g ab — N a N b is the intrinsic metric of X and D a is the metric covariant derivative of h. 
Thus the vanishing of K ab implies the vanishing of N a L ab h b , which gives the second condition. 

□ 

4 Conformal geodesies in Nil 

Nil is one of Thurston's list of eight geometries j4j. It admits a transitive 4-dimensional isometry 
group isomorphic to the 4-dimensional nilpotent Lie group. It can be considered with different 
signatures and since, as we shall see, the conformal geodesic equation can be reduced to quadra- 
tures, it gives information on the behaviour of conformal geodesies in the non-conformally-flat 
setting. 

We may take the metric to be 

g = dx 2 + edy 2 + (dz - \xdyf , (21) 
with e = ±1 in order to consider either Riemannian or Lorentzian forms. 

4.1 Riemannian Nil 

We'll start with Riemannian, so e = 1, and choose an orthonormal basis of one-forms: 

e 1 = dx, 6 2 = dy, 9 3 = dz- Xxdy, (22) 

with dual basis 

ei=d x , e 2 = d y + \xd z , e 3 = d z . (23) 
We calculate the connection one-forms as 
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and the nonzero Ricci tensor components as 

Rn = R22 = -R33 = -A 2 /2- 
We expand the unit velocity of a conformal geodesic in the chosen basis as 

u = aei + /3e 2 + 7e 3 , 

so that 

x = a, y = P, z- Xxy = 7, (24) 
and a 2 + /3 2 + 7 2 = 1. The acceleration is then calculated to be 

a = a\e\ + 0262 + 0363 

where 

ai = a + A/37 (25) 
02 = $ — Aa7 
a 3 = 7 

and the conformal geodesic equations then are 

di = --^702 - ^P a 3 - «|a| 2 - A 2 «7 2 (26) 

d 2 = ^701 + ^aa 3 - (3\a\ 2 - A 2 /?7 2 

03 = -^/3ai - -^aa 2 - 7|a| 2 + A 2 7(a 2 + /3 2 ) 

where |a| 2 = g a i,a a a h . From these, we claim that 

E ■= \a\ 2 - A 2 7 2 = constant. (27) 
This is just a matter of checking. To understand the constant E note that in general 

j t (g ab a a a b ) = 2L ab u a a b = j t (L ab u a u b ) - u a u b u c V a L bc , 

and for Nil the Ricci tensor satisfies the third Ledger Condition [6], equivalently the Killing 
tensor equation: 

V( a Rbc) = 0, 

so that \a\ 2 — L ab u a u b is constant since u a u b u c V a L bc = 0. E is a combination of this constant 
and a multiple of |u| 2 . 

Remark: The existence of the constant E provides a bound on |a| 2 , since necessarily 7 2 < 1. In 
turn this shows that, in Riemannian Nil, a conformal geodesic can be continued forever: there 
will not be any spirals in this case. 

Next note that the quantity 

J : = aiP - a 2 a + ^7 (28) 

is another constant of the motion. This is again a simple matter of checking, and this constant 
arises from the Killing vector d z analogously to the constants obtained in (|17p from Killing 
vectors. 
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It will follow from the existence of these two constants of the motion, that the equations are 
integrable. Note from the third of (|26p that 



7 = «3 = \{J ~ \>n) - l(E + A 2 7 2 ) + A 2 7 (l - 7 2 ), 



which integrates to give 



7 2 = -A 2 7 4 + (^A 2 - £) 7 2 + AJ 7 + Ci, 



in terms of a new constant C±, which we'll fix shortly. Suppose 7 known and parametrise a, (3 
by 

a = (l- 7 2 ) 1 / 2 cosx, /3 = (l- 7 2 ) 1/2 sin X . 
then the definition of J leads to 

X = (-J+^A 7 -A7 3 )(l-7 2 )^. (29) 

Now we express the acceleration in terms of x->l an d 7 to discover that C\ = E — J 2 . 
We have reduced the problem to quadratures: first for 7 we have 

7 2 = -A 2 7 4 + (^A 2 - E) 7 2 + AJ7 + E - J 2 =: F( 7 ); (30) 

then x is obtained from (|29p and the coordinates from (|24p . All quantities will be analytic in 
time, since the solution 7 of (|30|) is. Note in (j30l) that F(±l) = — (J=F 5 ) 2 so that, as desired, 7 
is confined within a range which can extend to 1 if J = A/2 or to —1 if J = —A/2, but not both; 
for general J neither limit is attained and for no orbit are both attained. 7 will be periodic in 
time but \ wm n °t m general, so that a, f3 also will not. In particular the orbits will not close 
in general. 

One example which can be solved explicitly is E = = J. Then 

7 2 = A 2 7 2 (jj-7 2 ), 

which is solved by 

7 = — secnu;, with uj = . 

' 2 2 

We may choose the origin of time so that u = ^^-t and solve for x to find 

3 sinh uj 

tan(x - Xoo) = , - 2 , 

1 — 2 sinn uj 

where Xoo is a constant of integration. Now as t — > 00 we have 7 — > and x tending to a 
constant. Thus a and /3 will tend to constants and the acceleration will tend to zero: this class 
of conformal geodesies tend to metric geodesies. It is easy to see that the metric geodesies which 
are also conformal geodesies fall into two classes: those with a = (3 = and 7 = ±1, and those 
with 7 = and a, (3 constants. It is the second set which are involved here. 

Remark: We can use the general solution to justify the claim made above, that the integral 
curves of a Killing vector which is not hypersurface orthgonal need not be conformal geodesies. 

To see this, consider the Killing vector 

K = /3d y + 7 d z 
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for constant /3,7 (which will shortly have the usual meanings). As a one-form this is 

K = (3dy + (7 - \/3x)(dz - Xxdy) 

so that K A dK = A(/3 2 — (7 — Xf3x) 2 )dx Ady Adz 7^ 0. Any trajectory of this Killing vector with 
x = has unit tangent u = K = (3e2 + 763, assuming that /3 2 + 7 2 = 1, and has acceleration 

a x = A/37, a 2 = a 3 = 0. 

Substituting into (j26|) we see that the trajectory is a conformal geodesic iff /3 7 = 0. So if we 
choose a Killing vector with ^7 7^ 0, then its integral curves are not conformal geodesies. 

4.2 Lorentzian Nil 

For the Lorentzian case, we take the metric to be 

g = dx 2 — dy 2 + (dz — Xxdy) 2 . 

We'll just consider space-like curves, so retaining (124ft a 2 — (3 2 + 7 2 = 1. Then (|26p becomes 

ai = d + A/37 (31) 
a 2 = /3 + Aery 
a 3 = 7 

and ([26]) becomes 

di = -^ 7 a 2 - - a(«i - a| + a i) + A 2 a7 2 (32) 

d 2 = ~^7 a i - ^aa 3 - P(a\ - a\ + a§) + X 2 ^ 2 

03 = ^/3ai - ^aa 2 - 7(0? - a\ + o§) - A 2 7(l - 7 2 ) 

The constants of the motion become 

E := af — a 2 + a 2 + X 2 j 2 = constant, (33) 

which does not now give a bound on the , and indeed the acceleration will not now be bounded 
in general, and 

J := a\/3 — 02« — -A7 = constant. 
The system reduces to quadratures again but now with 

7 2 = F( 7 ) := A 2 7 4 - (^A 2 + E) 7 2 + A J 7 + E + J 2 . (34) 

Note that F(±l) = ( J ± A/2) 2 , so 7 can go through ±1. 

The principal novelty in the Lorentzian case is that 7 can blow up in finite time: from (I34p 
once 7 is large we can have 7 ~ ±A7 2 so 7 ~ ±(A(io — *)) _1 for some to. 

To understand the conformal geodesies, we need to analyse the polynomial F, which is a 
quartic and may have 0, 2 or 4 zeroes. If it has none then 7 necessarily ranges from —00 to +00, 
or vice versa; if it has 2 then 7 comes in from ±00, meets a zero of F and goes back to ±00. 
If F has 4 zeroes then there are trajectories with unbounded 7 and others where 7 is confined 
between finite upper and lower limits, and so is periodic. 

For a and j3 we distinguish two cases: 
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• if 7 2 > 1, set a = (j 2 — l) 1//2 smhx, (3 = (j 2 — l) 1//2 coshx; 

• if 7 2 < 1, set a = (1 — 7 2 ) 1 / 2 coshx, /3 = (1 — 7 2 ) 1 / 2 sinhx- 
In either case, from the definition of J we obtain 

X = (J + ^A 7 -A 7 3 )(7 2 - 1)- 1 . 

This equation is singular at 7 2 = 1 but this is not really a problem. An unbounded 7 necessarily 
entails an unbounded /3, but it is not at this point clear how a behaves. A periodic 7 can give 
a nonzero secular change in x ° ver a period; if it does not then a, (3 are also periodic; otherwise 
they grow with one of a ± (3 — > 00 i.e. the velocity is asymptotically null. 
We again consider the special case E = = J. Then (|34|) becomes 

f = AV(7 2 " §), 

which integrates to give 7 = ^ sec u, where again oj = ^rpt with suitable choice of origin in t. 
Evidently 7 — > ±00 as t — > ±7r/(Av / 3): a finite time blow-up for the velocity. The integration 
for a and f3 leads to 

. _ (1-2 sin u) _ (1 + 2 sin u) 

a + /3 = Ci7— : f, a-/3 = C 2 - 



2(1 + sin w)' 2(1 -sine;)' 

where Ci and C 2 are fixed by data at t = 0, with C1C2 = 1. The combinations a ± /3 blow 
up at w = +vr/2 respectively. From (|24p it follows that the coordinate combinations x ± y will 
diverge in these limits too. It follows from (133f) that the components of acceleration and the 
scalar g e ja e af are also blowing up in finite time, which as we've seen cannot happen in the 
Riemannian Nil case. 

In any case in which 7 becomes large its clear from ()34j) that 7 runs away to infinity in 
finite time so that the terms in E and J in ()34f) become negligible and the case that we've just 
integrated can be seen to be typical. From the discussion of the zeroes of F we can say that, 
in any case in which 7 becomes large, the coordinates x and y diverge so there are no spirals, 
while if 7 is periodic then either a, f3 are also periodic or one of a ± f} diverges, whence so does 
one of x ± y and there are no spirals in this case. 

It remains to consider solutions with constant 7. If 7 = 70 = const, then the constants E 
and J and the normalisation of the velocity give 

2 a2 1 2 

or - P = 1 - 7 
a\a — 02/3 = 

ai/3 - a 2 a = J + -A70 

a\-a\ = E-X 2 j 2 

whence a\ = k{3,a,2 = ka for constant k = ( J + ^A7o)(7q — 1) _1 - Now from (I32p 

k/3 = a% = a + A/370, 

and similarly for /3, so that 

a = 0/3, /3 = £la, 

with Vt = k — A70. The solutions are exponentials (or constants if Vt = 0). Now the coordinates 
x, y will eventually be monotonic in time, diverging to ±00 and there are no spirals in this case. 
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An explicitly solvable special case with constant 7 is when a = j3; then 7 2 = 1 so we may 
suppose 7 = 1. The acceleration is 

o-i = 02 = ot + Aa, 03 = 
so that the acceleration is a null vector. The conformal geodesic equations collapse to 

a + Ad = — 7^A(d + Aa) + A 2 a, 

or 

3, 1 2 

a H — Aa A a = 0. 

2 2 

The solution is a sum of exponentials so that the velocity and acceleration are both unbounded 
(though bounded at any finite time). The orbit is not a spiral. 

In conclusion, in the case of the Riemannian or Lorentzian metrics on Nil, the conformal- 
geodesic equations can be reduced to quadratures and there are no spiralling conformal geodesies. 
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Appendix A: Conformal geodesies in the Berger Sphere 

We relegate to this Appendix the discussion of conformal geodesies in the metric known as the 
Berger sphere. This case also turns out to reduce to quadratures, but does not lead to new 
phenomena. The metric is 

where a 1 are the usual left-invariant one- forms on SU (2) , so that 

da 1 = a 2 Act 3 . 
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In Euler angles these one-forms can be taken to be 

a 1 + ia 2 = e^(d6 - i sin 9d(f)), a 3 = dip + cos Qdcfi, 
when the dual basis of vector-fields is 

e # (<%--^-^ + icot^), d^. 

Clearly the metric becomes the round metric on S 3 when A = 1. For the curvature calculation, 
we take the orthonormal basis of forms to be 

Q^ = a\ 6 2 = o 2 1 6 3 = \a\ 

then the connection one-forms are 

u\ = C9 3 , w 2 3 = A9\ uj 3 1 = AO 2 , 
where A = A/2 and C = (2 — A 2 )/2A. The nonzero Ricci components are 

Ru = R22 = 1 — A 2 /2, i?33 = A 2 /2, 

so the curvature scalar is R = 2 — A 2 /2. 

We expand the velocity u and acceleration a in terms of the basis (ei, 62,63) of vector fields 
dual to the basis (9 l ) as 

u = aei + /3e 2 + 7e 3 , 
a = a\e\ + 0262 + 0363, 
then the acceleration has components 

ax = a + (C - A)/3j 
a 2 = $-(C- A)^a 
a3 = 7, 

and the conformal geodesic equations become 

a'i = —6*702 + A(3(i3 — a(a 2 + a 2 4- a 2 ) + i?na(l — a 2 — (3 2 ) — R^aj 2 
a 2 = -Aaa 3 + C^m - /3(a 2 + a\ + a§) + i?n/3(l - a 2 - /3 2 ) - i? 3 3/37 2 
a 3 = - A/3ai + Aaa 2 - 7(0? + a 2 + 03) " ^n7(« 2 + /? 2 ) + ^337(1 - 7 2 ) ■ 

As before we seek constants of the motion, and find 

J := a\(3 — (12OL — A'j = constant, 

which is analogous to (|28p and obtained from the Killing vector d^, and 

E := \a\ 2 + 7 2 (i?n - #33) = a\ + a% + a\ + (1 - A 2 )7 2 = constant, 

which is analogous to (|27|) . E serves to bound the magnitude of the acceleration, and exists 
because the Ricci tensor is again a Killing tensor. We derive an equation for 7 alone: 

7 = (1 - A 2 )(2 7 3 - 7) - (E + yl 2 ) 7 - AJ 

We shall integrate this shortly, but before that set 

a = Rcosx, P = Rsinx, where R = (1 — 7 2 ) 1 / 2 , 
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for then the definition of J leads to 

. J + A 7 -(C-^) 7 (l-7 2 ) 
X ~ 1 - 7 2 

and then 

( 7 ) 2 = F( 7 ) := (1 - A 2 )( 7 4 - 7 2 ) - (£? + .4 2 ) 7 2 - 2AJ^ + E-J 2 . 

Note that = — (J±^4) 2 , so that the evolution forces 7 2 < 1. The general discussion is 

now just like the case of Riemannian Nil: since the acceleration is bounded, the orbits can be 
extended indefinitely but will not close in general. 

The special solutions with E = = J can be found explicitly: 

( 7 ) 2 = (1 - A 2 )( 7 4 - 7 2 ) - AW = 7 2 ((1 - A 2 ) 7 2 - (1 - ^A 2 )). 

These do not exist for 4/3 > A 2 but for A 2 > 4/3 we find 

7 = /c sech/ui, with k 2 = — — — , n = k(\ 2 - 1) 1/2 . 

A — 1 

Now 7 , 7 — > as t — > ±oo, so that in these limits the acceleration goes to zero and the conformal 
geodesic asymptotes to a metric geodesic. The metric geodesies which are conformal geodesies 
form two classes: those with 7 2 = 1 so that a = = /3; and those with 7 = and a, (3 constants. 
It is the second family which are serving as asymptotes here. Brief calculation shows this class 
are closed (as are the other class) but this is not true for all geodesies. 

As with Nil there is no spiralling, which may be a reflection of the integrability. 

Appendix B: Conformal geodesies in the Schwarzschild solution 

The Schwarzschild solution is a familiar Lorentzian Einstein metric, so that the conformal 
geodesic equations imply that the magnitude of the acceleration is constant, but there is no 
reason to expect the components of the acceleration to be bounded. This case is worth studying 
as there are two constants of the motion arising from the Killing spinor, although this is not 
sufficient to reduce the system to quadratures. 

The metric in the region outside the horizon can be taken to be 

- V 2 dt 2 + \dr 2 + r 2 {d6 2 + sin 2 Odd) 2 ), (Bl) 

with V 2 = 1 — 2m/r, (so we are restricting to r > 2m). We take the orthonormal basis of 
one-forms as 

6° = Vdt, 6 1 = V~ 1 dr, 6 2 = rd9, 6 z = rsmed(j). 
Then expand in this basis: 

u = a9° + /39 1 + 7 # 2 + 56 3 
a = a 6° + aiO 1 + a 2 2 + a 3 6 3 
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and calculate 



a = ao — aj3V' 

= a x - aV + (7 2 + 5 2 )V/r 

7 = 02 — fl^/V/r + 5 2 cot 6/r 

6 = 03 — j38V/r — j5 cot 6/r 

clq = —eAa — aaiV' 

a\ = —eA[3 — aaoV + ( , ya2 + Sas)V/r 

a2 = — eAj — ja%V/r + ^03 cot 6/r 

03 = — eA5 — 5a\V/r — 5a2 cot 6/r 

where 

A = — ag + a\ + a\ + 03. 

This is the norm of the acceleration which is constant but may have either sign, and e is -1 for 
time-like u and +1 for space-like, so that 

e = -a 2 + /3 2 + 7 2 + 5 2 , 

and of course 

-aa + /5a\ + 702 + Sa 3 = 0. 

For time-like u, the acceleration is space-like and A must be positive. The following can be 
checked to be constants of the motion: 



r(/3ao — aai) — eaV, r{5a2 — 703)- 
These are the real and imaginary parts of the quantity 

Q = u a a b F ah - eK h u\ (B2) 
where = ojab^a'B' ls obtained from the Killing spinor loab which satisfies 

V ' AA'^BC = ^K A ^ B e C )A 

with K a = Kaa 1 the time-like Killing vector, unit at infinity. There do not seem to be constants 
of the motion related directly to the (four linearly-independent) Killing vectors of the solution 
so it is interesting that two arise from the Killing spinor. From the definition, it is simple to 
check that Q in (lB2j) is a constant of the motion, and will be as well for the Kerr solution, which 
likewise possesses a Killing spinor. 

We may explore across the horizon by introducing A = V(a + j3)/2,v = (a — fi)/V and 
correspondingly ai = V(ao + ai)/2, a^v = (ao — a i)/V ■, which follow by introducing the two null 
vectors L = (eo + e-xj/V, N = V(eo — e\)/2 into the basis in place of eo, e\. The corresponding 
null coordinates are obtained by solving 

dr dr 
du = dt-—K, dv = dt + -pr, (B3) 
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and the system regular at the horizon is 



V 2 

A = a L + — (2\u + e) - VV'Xv 
2r 

v = a N --(2\v + e) + VV'v 2 

r 

7 = a 2 --(\--vV 2 ) + 5 2 cot9/r 

IT 

S 1 

5 = o 3 - -(A- -vV 2 ) - 7 (5 cot# jr 

V _■ 

at = — eAA — VV v(il + — — {vcll + Aajv) 

oat = — eAi/ + VV'vcin (vai + Aajy) 

r 

7 1 2 

d2 = —eA'y (at — -a^V ) + 5a^cot6/r 

IT — 

S 1 

0,3 = — eA<5 (aL — -^o^nV 2 ) — ^2 cot 9/r 

IT — 



The first constant of the motion can be written 

r(\a N - va L ) - e(A + l;V 2 v)- 

Without integrability we are restricted to integrating some special cases. The case of radial 
time-like conformal geodesies, in the system which extends across the horizon, simplifies to 

A = 2^ 2 ^ = 2 gX ' v = ^ = e ~ X ^ ° L = 2 agX ' aN = ~ ae ~ X ^ 



with constant a, A = a 2 and 
which integrates to give 



X = a - VV'e- x , 



e x + V 2 e~ x = 2ar + C 



with constant of integration C. If a 7^ it is convenient to set C = —2aro when we can derive 
the following equation for r: 

r 2 = -V 2 + a 2 (r - r ) 2 = a 2 (r - r ) 2 - 1 + ^, 

from which we obtain information about turning points. Note that r is zero where V 2 = 
a 2 (r — ro) 2 which can be solved graphically. If vq > 2m then there are two solutions to this 
equation, both greater than 2m; for one the conformal geodesic comes in from past null infinity 
Z~, has a turning point and goes out to future null infinity I + ; for the other the conformal 
geodesic comes from the past singularity in the Penrose diagram, crosses r = 2m and reaches a 
maximum value of r, before recrossing r = 2m and ending at the future singularity. Reducing 
ro, these solutions will persist until a critical value of ro is reached, when the graphs of V 2 and 
a 2 (r — ro) 2 are tangent; now there is a trajectory of the Killing vector d/dt which is a conformal 
geodesic (as is to be expected from Theorem 3.1), and there are incoming and outgoing radial 
conformal geodesies asymptotic to it. Reducing ro still further, r has no zero and we obtain 
incoming or outgoing radial conformal geodesies with no turning point. 
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Another integrable case is that of solutions confined to the equatorial plane. Set a = 7 = 
o-q = a.2 = cos 9 = 0. There are circles like this with (3 = too which are trajectories of the 
d/dcf> Killing vector, but more generally we can set 



P = cos x, ^ = sin x, 01 = -q sin x, a 3 = q cos x, 
then q = x + (Y/ r ) sinx = constant. This implies 

(rsinx)' = qr cos x = qrr/V, 

which is solvable as 

r 2 = V 2 {1 - {G + Cf/r 2 ) 

with constant C and G' = qr/V. The integral could be done explicitly, and there will be 
solutions which do not wind around the centre as well as ones which do, but the curves will not 
close or spiral. 
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